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ABSTRACT 


Analytical  and  experimental  work  on  microsegregation  and 
dendrite  structure  in  a  series  of  iron  base  alloys  is  summarized. 
Experimental  work  is  on  iron-nickel  alloys,  iron-4 %  phosphorous 
alloy,  iron-25%  copper  alloy,  and  iron-nickel  carbon  alloys. 

A  new  computer  analysis  for  microsegregation  is  presented, 
with  detailed  numerical  results  given  for  iron-26%  nickel  alloy, 
iron-10%  nickel  alloy,  and  iron-4%  phosphorous  alloy.  Results 
agree  reasonably  well  with  experiment  provided  the  dendrite  arm 
spacing  used  in  calculation  is  less  by  a  constant  "correction 
factor",  g,  than  that  measured  (for  example  g  ~  0.40  for  cylinder 
model,  secondary  dendrite  arm  spacing).  It  is  concluded  the  need 
for  the  correction  factor  arises  primarily  because  the  model 
employed  for  calculation  does  not  account  for  "ripening"  during 
solidification . 

Structures  of  the  various  alloys  studied  are  presented;  in 
some,  the  dendrite  morphology  is  strongly  plate-like;  in  others, 
it  is  "fibrous".  Over  a  range  of  cooling  rates  from  10-3oC/sec 
to  nearly  10^°C/sec  dendrite  arm  spacing  of  iron-26%  nickel 
alloy  is  linearly  related  to  a  single  power  of  cooling  rate.  It 
is  concluded  that  the  dependency  results  from  the  predominant 
effect  of  "coarsening"  in  determining  final  dendrite  arm  spacing. 
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CHAPTER  I 


INTRODUCTION 


Research  has  been  conducted  at  Massachusetts  Institute  of 
Technology  on  microsegregation  in  iron-base  alloys  since  late 

1-4 

1962.  Results  have  been  summarized  in  annual  reports  and 

portions  have  been  abstracted  and  published  as  technical 
5-7 

papers  ;  a  series  of  technical  papers  are  now  being  prepared 
covering  the  later  work. 

General  aim  of  the  program  has  been  to  obtain  a  sound 
engineering  understanding  of  factors  influencing  the  form  and 
distribution  of  microsegregation  in  cast  steel.  Specific 
aspects  of  the  work  have  included: 

1.  Study  of  dendrite  morphology  and  dendrite  arm  spacing 
in  iron-base  alloys. 

2.  Analytical  and  experimental  studies  on  solute 
distribution  (microsegregation)  in  binary  iron-base 
alloys . 

3.  Analytical  and  experimental  studies  on  solute 
distribution  (microsegregation)  in  ternary  iron-base 
alloys . 

4.  Analytical  and  experimental  studies  on  homogenization 
kinetics  in  iron-base  alloys. 

This  report  summarizes  the  last  part  of  research  in  this 
program  devoted  primarily  to  study  of  solidification  of  the 
primary  (dendrite)  phase  in  iron-base  alloys.  Current  and 
planned  research  deals  primarily  with  second  phases.  Present 
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studies  emphasize  the  formation,  distribution,  and  morphology 
of  second  phases  (especially  non-metallic  inclusions)  precipi¬ 
tating  during  solidification. 

Chapter  2  of  this  report  summarizes  results  of  computer 
studies  on  microsegregation  in  two  different  binary  alloy 
systems.  Iron-nickel  was  chosen  as  an  example  of  a  system  in 
which  essentially  complete  solid  solubility  is  obtained,  and 
iron-phosphorous  as  an  example  of  a  system  containing  a  eutectic. 
Portions  of  the  work  on  iron-nickel  were  presented  earlier  ,  but 
the  bulk  of  this  chapter  deals  with  results  of  a  new  computer 
analysis  of  microsegregation. 

Chapter  3  summarizes  briefly  the  melting  and  casting 

procedures  used  in  experimental  portions  of  this  work  and 

Chapter  4  presents  results  of  studies  on  dendrite  morphology 

and  dendrite  arm  spacing  measurements.  Again,  portions  of 

earlier  work  on  iron-nickel  are  presented  in  condensed  form  for 

continuity,  but  most  of  the  work  is  new  including  the  important 

correlation  of  Figure  4-7  which  shows  that  dendrite  arm  spacing 

is  linearly  related  to  a  single  power  of  cooling  rate  in  iron- 

25%  nickel  alloy  over  at  least  a  range  of  cooling  rates  from 
“3  6  ' 

10  °C/sec  to  10  °C/sec.  The  former  rate  was  obtained  in  a 
controlled  cooling  furnace  and  the  latter  in  a  "splat  cooling" 
apparatus.  Intermediate  cooling  rates  were  obtained  in 

i 

unidirectionally  solidified  ingots.  Reasons  for  the  dependency 
of  dendrite  arm  spacing  on  cooling  rate  are  discussed. 
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Chapter  5  summarizes  measurements  of  microsegregation 
made  in  binary  iron-base  alloys  over  the  last  several  years 
and  compares  these  measurements  with  the  predictions  of  the 
new  computer  analysis  of  Chapter  2. 
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CHAPTER  2  -  ANALYSES  OF  MICROSEGREGATION 
IN  BINARY  IRON-BASE  ALLOYS 

A.  Summary  of  Analyses  Employed 

In  this  section,  a  summary  is  given  of  three  methods 
employed  for  calculation  of  microsegregation  in  binary  iron- 
base  alloys.  The  first  two  of  these  were  developed  in  an 
earlier  research  period  and  reported  elsewhere1'2 3.  A  new 
method  of  analysis  was  developed  in  this  work  with  the  aims 
of  (1)  treating  extensive  diffusion  in  the  solid  with  greater 
accuracy,  and  (2)  treating  the  case  of  cylindrical  dendrite 
geometry . 


Computations  are  for  a  local  characteristic  volume  element. 
Assumptions  of  all  methods  of  analysis  employed  include  the 
following:1 

1.  There  is  negligible  undercooling  before  nucleation,  or 
from  curvature  or  kinetic  effects. 

2 .  There  is  no  mass  flow  in  or  out  of  the  volume  element 
considered.  Such  flow  might  occur,  for  example,  from 
movement  of  liquid  by  convection  or  to  feed  shrinkage, 
from  movement  of  solid  during  solidification,  or  from 
diffusion . 

3.  Diffusion  in  the  liquid  within  the  volume  element 
considered  is  complete. 

The  equilibrium  partition  ratio,  k,  applies  at  the 
interface . 


4. 
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The  three  methods  of  calculating  microsegregation  are 
summarized  below  and  their  differences  highlighted. 

Method  I :  Analytic  Solution 


For  the  solidification  model  described  by  the  above 

conditions,  the  redistribution  of  solute  during  solidification 

‘ 

initially  was  described  analytically  by  making  the  following 

I 

assumptions : 

1.  The  dendrite  model  is  platelike. 

2.  The  rate  of  advance  of  the  interface  is  either  linear 
or  parabolic. 

dAi 

Linear:  |  A^«t  and  ^  =  constant  (2-1)* 

i 

Parabolic:  A.  «/t  and  -rr- —  =  constant  (2-2) 

|  ’  x  dt  /t 

3.  Diffusion  in  the  solid  is  sufficiently  small  that  the 
diffusion  does  not  affect  significantly  the  composi¬ 
tion  in  the  solid  at  the  interface 


dC* 

s 

d  A. 

r 


(2-3) 


Now  the i differential  materials  balance 

|  d(fLCL)  +  d(fsCs) .=  0  (2-4) 

i 

| 

is  evaluated  for  the  above  assumptions  using  the 
requirement  that  the  total  amount  of  solute  entering 
the  solid  due  to  diffusion  in  the  solid  is  equal  to 
the  diffusion  flux  at  the  interface.  Then 


*  Symbols  used  jin  the  text  are  defined  in  Appendix  A. 
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A  A. 
1 

dT” 


constant 


(CT -C*)  df 
L  s  s 


ctdC*+  (1-f  )  ac 


constant 

/§■ 


(CL--Cs)dfS 


where 


a 


D  t. 
s  f 


2afsdC*+(l-fs)ACL 


(2-5) 

(2-6) 

(2-7) 


For  constant  partition  ratio,  k,  and  constant  solid 
diffusion  coefficient,  D^,  the  differential  equations 
(2-5)  and  (2-6)  are  easily  integrated  to  obtain  simple 
analytic  solutions.  For  the  more  general  case  of 
variable  k  and  Dg ,  the  differential  equations  (2-5)  and 
(2-6)  are  solved  numerically.  During  the  isothermal 
solidification  that  may  occur  at  the  eutectic  temperature, 
solid  diffusion  diminishes  the  amount  of  material  of 
eutectic  composition.  Method  I  may  be  used  to  commute 
the  change  at  the  eutectic  temperature  in  the  eutectic 
fraction.  Method  I,  however,  has  not  been  applied  to  the 
solute  redistribution  that  occurs  during  cooling  after 
solidification . 


Method  I  overestimates  the  amount  of  solid  diffusion  and, 
thus,  underestimates  the  extent  of  microsegregation.  Errors 
become  serious  when  ak  approaches  1.  Finally  Method  I  predicts 
solid  composition  at  the  solid  liquid  interface,  C*,  but  does 
not  yield  a  description  of  the  solute  distribution  within  the 


dendrite . 
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Method  II:  Numerical  Analysis 

To  account  for  solid  diffusion  after  solidification  and  to 
achieve  a  complete  description  of  the  solute  distribution  within 
the  dendrite  both: during  and  after  solidification,  a  simple 
numerical  method  was  employed.  For  Method  II 

(1)  The  dendrites  again  are  considered  platelike 

(2)  The  rate  jof  advance  of  the  interface  is  either  linear 
or  parabolic  (equations  2-5  and  2-6) 

(3)  The  solid  composition  at  the  liquid  solid  interface, 

C*,  at  a  particular  fraction  solid,  f  ,  is  obtained 

*  S 

by  Method  I . 

i 

To  evaluate  splid  diffusion  during  solidification  Fick's 

second  law  of  diffusion  is  transformed  by  finite  differences  to 

the  approximate  expression 

Cc<  +  C  (M  -  2)  +  C 

C  =  n+r,k-i  t,k-i _ .  ,, 

S  .  ,  M  '  **  '  -  > 

.  M  A  A2  ! 

where  M  =  i  (2-9) 

I 

i 

This  equation  is  evaluated  iteratively  over  distance  (from  A  =  0 
to  A^)  and  over  time  (from  t  =  0  to  tf)  for  the  initial  and 
boundary  conditions 


when  t  =  0,  A.  =  0  !  and  C  (A.)  =  kC 
1  s  i  o 

9C  j 

at  A  =  0 ,  -rr—-  =  0  and  at  A  =  A.  ,  C  (  A.  )  =  C* 

9  A  i  i  '  s  '  l  s 

I  Ai 

where  C*  is  given  by  the  Method  I  for  f  =  . 

S  S  jL 


(2-10) 
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Method  II,  also,  may  be  used  to  calculate  the  changes  in 
microsegregation  at  the  eutectic  temperature.  In  addition. 
Method  II  is  readily  applied  to  the  computation  of  the  solute 
redistribution  on  cooling  after  solidification  by  iteratively 
evaluating  equation  (2-8)  using  as  the  initial  condition  the 
solute  distribution  at  the  actual  solidus  temperature  and  with 
the  boundary  conditions  for  a  eutectic  alloy; 


at  X.  =0,  =  0;  and  at  A  =  A.  , 

1  d  A  X 

and  for  a  single  phase  alloy; 

9Cs 

at  A.  =0,  =  0;  and  at  A  =  A.  , 

i  3  A  l 


C  (A.) 
s  i 


=  Cg (max.  solubility) 
.  (2-11) 

0  (2-12) 


Method  II  is  more  accurate  than  Method  I,  yet  because  it 
depends  on  the  results  of  Method  I  for  interface  composition, 
it  to  cannot  be  used  when  there  is  a  large  amount  of  solid 
diffusion.  Method  II  underestimates  the  amount  of  solid  diffu¬ 
sion  and  solute  mass  is  not  conserved.  Use  of  a  correction 

2 

factor  described  previously  ,  fortuitously,  allows  the  accurate 
use  of  Method  II  for  large  amounts  of  solid  diffusion.  However, 
the  modified  method  is  not  straightforward  and  is  time 
consuming . 


Method  III:  Mass  Balance  and  Finite  Difference 

Method  III  is  the  new  numerical  technique  and  computer 
program  that  has  been  developed  to  account  readily  for  large 
amounts  of  solid  diffusion  and  to  allow  for  a  variation  of 
dendrite  geometry.  In  this  method: 


(1)  The  dendrite  model  may  be  taken  as  platelike, 
cylindrical,  or  spherical. 

(2)  The  rate  of  advance  of  the  interface  may  be  linear, 
parabolic  or  governed  by  thermal  data, 

(3)  The  composition  in  the  solid  at  the  interface  is  set 
equal  to  kC^  where  is  determined  by  mass  balances 

i 

after  every  interval  of  solidification. 

Briefly,  diffusion  within  the  dendrite  is  computed  using 
the  finite  difference  technique,  i.e.  by  iteratively  evaluating 
the  algorithm  ( 2  —  8 )  for  0  <  A  <  A^ .  Then  the  mass  balance, 
equation  (2-4) ,  is  used  to  compute  the  new  composition  of  the 
liquid.  An  amount  A  A  of  solid  is  considered  to  form  at  the 
composition  that  is  in  equilibrium  with  the  liquid,  kCT ,  and 
the  process  is  repeated.  This  method  may  also  be  used  to 
compute  solute  redistribution  at  the  eutectic  transformation  and 
on  cooling  after  solidification.  The  details  of  the  numerical 
method  are  given  in  Appendix  B. 

In  summary.  Method  III  is  a  straightforward  technique 
readily  used  to  describe  solute  distribution  within  a  dendrite 
during  solidification  and  on  cooling  after  solidification.  The 
method  is  accurate  even  for  large  amounts  of  solid  diffusion 
and  may  be  used  for  cylindrical  as  well  as  platelike  dendrite 
geometries . 
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B .  Alloy  Systems  and  Diffusion  Data  Used  as  Examples 

Two  binary  alloy  systems  were  used  in  this  study  as 
examples  for  calculation.  The  first  of  these  is  the  iron- 
nickel  system.  In  this  system  alloys  above  about  6  per  cent 
nickel  solidify  completely  as  a  single  phase,  austenite.  In 
all  iron-nickel  alloys,  the  composition  of  interdendritic  spaces 
is  expected  to  approach,  as  a  limit,  the  composition  of  the 
liquidus  minimum  at  68%  nickel.  Figure  2-1.  The  second  alloy 
system  used  as  example  is  the  iron-phosphorous  system,  which 
has  a  simple  eutectic  at  the  iron-rich  end  with  limited  solid 
solubility,  Figure  2-2. 

Data  for  the  diffusion  coefficient  used  in  calculations  for 

.  3 

iron-nickel  alloys  are  from  Goldstein  who  showed  the  diffusion 
coefficient  for  the  y  phase  varies  with  both  temperature  and 
nickel  content  according  to  the  relation,* 

38,380  -  5.85  X  . 

Dg  =  exp  [0.0519  XNi  +  1.51]  exp  [ - ^ - — ]  (2-13) 

where  XrT.  —  atomic  per  cent  nickel 
Ni 

T  =  temperature,  degrees,  Kelvin 

Data  for  the  diffusion  coefficient  used  in  the  calculations 

4 

for  the  iron-phosphorus  system  were  taken  from  Seibel  who 

*  Goldstein  reports  the  above  expression  to  be  in  agreement  with 
his  measured  values  of  Ds  within  ten  per  cent  for  nickel  con¬ 
tents  up  to  50  per  cent  in  the  temperature  range  of  1000°C  to 
1288°C.  In  the  work  equation  (2-1)  is  employed  for  composi¬ 
tions  up  to  68  per  cent  nickel  and  in  the  temperature  ranges 
800°C  to  1469. 5°C. 
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lists  the  following  expression  for  phosphorus  diffusion 
in  a : 

n  —o  an  ,,  r  55000-, 

Dg  =  2.90  exp  [ - - - ]  (2-14) 

where  R  =  gas  constant. 

In  this  work/  solutions  are  presented  only  for  a  linear 
volumetric  growth  rate  and  for  platelike  and  cylindrical 
dendrite  geometries. 

C.  Results  of  Calculations ,  Iron-26%  Nickel  Alloy 
Method  I : 

I 

Details  of  this  analysis  (the  simplest  employed  herein) 

have  been  given  earlier.  Results  for  iron-26%  nickel  alloy 

are  given  in  Figure  2-3,  which  plots  solid  composition  at  the 

interface,  C*,  versus  fraction  solid,  f  ,  for  various  values 
s  s 

of  ,  where: 


At  n  -  0,  diffusion  in  the  solid  is  negligible  and  the 
calculated  curve  is  that  resulting  from  the  Scheil  equation 
(with  variable  k)  .  ;  At  n  =  diffusion  in  the  solid  is  complete 
and  the  result  is  that  obtained  from  the  equilibrium  lever  rule. 

i 

! 

i 

If  diffusion  in  the  solid  is  slight  (n  <<  1) ,  then  the 
curves  of  Figure  2-3  also  represent  approximately  the  solid 
concentration,  Cg ,  that  would  be  measured  from  a  microprobe 
trace  across  the  pfatelike  dendrite  arms.  This  trace  would 
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start  at  the  center  of  the  arm  (f  =  0)  and  end  at  the  center 

O 

plane  between  the  two  dendrite  arms  (f  =  1) .  It  is  clear, 
however,  that  Figure  2-3  cannot  be  used  in  this  way  as  n  -  ro, 
for  example,  the  interface  composition  C*  is  as  shown  in  Figure 
2-3.  However,  the  final  solid  composition,  C  ,  would  be  26% 
nickel  at  all  values  of  f  . 

Method  II : 

This  finite  difference  method  (Method  II,  Appendix  A)  was 
employed  previously  to  give  a  better  description  of  composition 
across  dendrite  arms.  For  substantial  diffusion  such  as  is 
clearly  present  in  iron-26%  nickel,  this  method  has  been  super- 
ceded  by  Method  III,  the  "mass  balance  technique".  Some  results 
of  this  technique  are,  however,  discussed  below  for  comparison 
with  Method  III. 

Method  III: 

Typical  results  of  Method  III  are  given  in  Figure  2-4. 

Figure  2-4  shows  solid  composition,  C  ,  versus  fractional 
distance,  A/L,  along  a  dendrite  arm  (after  solidification)  for 
a  platelike  dendrite  model,  assuming  r\  =  4  x  10^  and  linear 
advance  of  the  solid-liquid  interface  with  time.  Also  shown  in 
Figure  2-4  are  results  for  the  cylindrical  dendrite  model. 
Segregation  in  the  cylindrical  model  is  significantly  less  than 
that  in  the  plate  model.  Comparison  of  Figure  2-4  with  2-3 
shows  that  in  the  plate  model,  n  =  4  x  10^,  results  of  Methods 
I  and  III  are  closely  similar.  The  results  of  Method  II  and 
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Method  III  are  compared  in  Figure  2-5  which  gives  "Segregation 
Ratios"*  at  the  solidus  temperature.  S'  versus  k  calculated  by 
Methods  II  and  III  (plate  model) .  Segregation  ratios  at  all 
values  of  n  are  closely  similar  (but  slightly  less  for  Method 
III) *. 

Figure  2-6  shbws  local  maximum  and  minimum  solute  contents 

at  the  actual  solidus  temperature,  and  CV,  plotted  individually 

versus  n.  Also  shown  are  the  local  maximum  and  minimum  solute 

contents  at  room  temperature,  C^,  C^.  Calculation  of  room 

temperature  values  was  based  on  cooling  curve  data  obtained  from 

unidirectional ly  solidified  ingots  as  discussed  in  Chapter  3; 

and  the  method  of  calculation  is  described  in  Appendix  B.  It  is 

evident  in  Figure  2-6  that  the  ratio  of  CfVC0  is  substantially 

Mm  x 

less  than  C^/C^  predicted  as  a  result  of  diffusion  in  the  solid 
during  cooling. 

Figure  2-7  shows  the  variation  of  both  S'  and  S°  (where 
S°  =  C^/C^)  with  n  for  plate  and  cylinder  dendrite  models. 

*  Segregation  Ratio,  S,  is  defined  as: 


CM  =  maximum  composition  of  the  solid  within  regions  the 
order  of  a  dendrite  arm  spacing 

=  minimum  composition  of  the  solid  within  regions  the 
order  of  a  dendrite  arm  spacing 

The  quantities  are  primed  (S',  C'  C')  when  taken  at  the  non¬ 
equilibrium  solidus  temperature  and  written  as  S°,  C° ,  C° 
when  taken  at  room  temperature.  M  m 


Lowest  segregation  is  for  the  cylinder  model  at  room 
temperature . 
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D .  Results  of  Calculations  -  Iron-10%  Nickel  Alloy 

Figures  2-8  and  2-9  show  summary  results  of  calculations 
on  iron-10%  nickel  performed  similar  to  the  foregoing  (Method 
III) .  The  following  general  observations  may  be  made  from 
these  figures: 

(1)  Significant  homogenization  occurs  in  this  alloy  after 
solidification  and  during  cooling  to  room  temperature 
for  all  values  of  n  plotted  (as  in  iron-26%  nickel 
alloy) . 

(2)  The  cylindrical  dendrite  model  predicts  less  segregation 
at  all  values  of  n  than  does  the  platelike  model  (as  in 
iron-26%  nickel  alloy) . 

(3)  Calculated  values  of  S'  and  S°  are  significantly  higher 
for  iron-10%  nickel  alloy  than  for  iron-26%  nickel 
alloy. 

E,  Results  of  Calculations  -  Iron-4%  Phosphorus  Alloy 

Figures  2-10  and  2-11  show  results  of  calculations  on  the 
iron-phosphorous  system  (see  phase  diagram,  Figure  2-2) .  This 
alloy  system  contains  a  eutectic  at  10.2%  P,  and  the  limit  of 
solid  solubility  at  the  eutectic  temperature  is  2.8%. 

Figure  2-10  shows  calculated  weight  per  cent  eutectic  versus 
n  for  iron-4. 0%  phosphorous  alloy  at  the  solidus  (eutectic) 
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temperature,  and  for  room  temperature.  Note  the  extent  of  solid 
state  diffusion  after  solidification  is  much  less  in  this  alloy 
than  in  the  iron-nickel  alloys,  primarily  because  the  concentra¬ 
tion  gradients  in  the  solid  are  much  less. 

Figure  2-11  shows  how  predicted  minimum  composition  (at  the 
center  of  dendrite  arms)  should  vary  with  n.  This  minimum  goes 
from  kCQ  at  low  values  of  n  to  the  maximum  solid  solubility  (2.8%) 
at  higher  values  bf  ri , 

F ,  Conclusions 

1.  A  new  computer  program  was  written  to  simulate 
solidification  and  post-solidification  heat  treatments  by  using 
numerical  analysis  techniques.  The  program  may  be  applied  to 
eutectic  alloys  as  well  as  those  in  which  (a)  the  final  solidi¬ 
fication  temperature  and  composition  are  variable;  and 

(b)  substantial  diffusion  occurs  in  the  solid  both  during  and 
after  solidification.  Characteristics  of  the  program  are  that 
(a)  it  precisely  conserves  solute  mass,  (b)  it  permits  assump¬ 
tion  of  either  a  cylindrical  or  platelike  dendrite,  (c)  growth 
rate  in  the  solidification  range  may  be  linear,  parabolic,  or 
determined  by  a  cooling  curve  or  cooling  expression,  (d)  any 
desired  cooling  curve  or  cooling  expression  may  be  employed 

i 

from  the  solidus  temperature  to  room  temperature. 

2.  For  values  of  r\  (local  solidification  time  divided  by 
the  square  of  half!  the  dendrite  arm  spacing)  that  will  later 
be  seen  to  be  experimentally  obtained  for  iron-nickel  alloys: 
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(a)  Significant  diffusion  occurs  in  the  solid  during 
solidification . 

(bj  Significant  diffusion  occurs  in  the  solid  after 
solidification  and  during  cooling  to  room 
temperature. 

3,  For  a  given  value  of  n ,  predicted  segregation  ratios, 
S'  (at  the  solidus  temperature)  and  S°  (at  room  temperature) 
are  (a)  less  for  iron-26%  nickel  than  for  iron-10%  nickel,  and 
(b)  less  in  both  alloys  for  cylindrical  model  than  for  the 
plate  model. 

4.  In  iron-phosphorous  alloys,  residual  eutectic  and 
minimum  concentration  at  the  centers  of  dendrite  arms  depend 
on  n .  At  a  given  value  of  g ,  extent  of  diffusion  in  the  solid 
is  much  less  than  for  the  iron-nickel  alloys  studied. 
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Fe-Ni  PHASE  DIAGRAM 


WEIGHT  PERCENT  NICKEL 


TEMPERATURE  (°C) 


FRACTIONAL  DISTANCE  ALONG  DENDRITE 

Figure  2-4.  Predicted  solute  distribution  at  the  solidus  temperature 
for  an  iron-26  per  cent  nickel  alloy  with  r  =  4  x  10^ 
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Figure  2-5.  Predicted  segregation  ratios  at  the  solidus,  S  ,  versus 
for  iron-26  per  cent  nickel  alloy,  Plate  model. 


Figure  2—6.  Predicted  maximal  and  minimum  solute  contents  as  a 

function  of  ri  for  an  iron-26  per  cent  nickel  alloy  at 
the  solidus  and  at  room  temperature.  Plate  model. 
Method  III. 


Fe-10%  Ni 


Figure  2-8.  Predicted  maximum  and  minimum  solute  contents  as  a  function 
of  r\  for  iron-10  per  cent  nickel  alloy  at  the  solidus  and 
at  room  temperature.  Plate  model.  Method  III* 


Fe  -26% Ni 


Figure  2-7.  Predicted  values  of  the  segregation  ratio  at  the  solidus 
S',  and  at  room  temperature,  S°,  as  a  function  of  n  for 
iron-26  per  cent  nickel  alloy.  Method  HI. 


Fe  -  4.0  %  P 


Figure  2-10.  Weight  per  cent  eutectic  versus  n  for  an  iron-4.0  per  cent 
phosphorus  alloy.  Plate  Model.  Method  III. 
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CHAPTER  3  -  CASTING  PREPARATION  AND  THERMAL  ANALYSIS 

i 

A.  Unidirectionally  Solidified  Plate  Castings 

A  number  of  alloys  were  prepared  in  the  iron-nickel  iron- 

I 

phosphorus,  and  iron-copper  binary  systems,  and  the  iron- 
nickel-carbon  ternary  system  for  study  of  microsegregation  and 

dendrite  morphology.  A  list  of  the  alloys  studied,  the  charge 

| 

compositions,  and  the  results  of  chemical  analyses  are  given 
in  Table  3-1.  i 

All  heats  were  induction  melted  in  clay  graphite  crucibles 
with  rammed  linings  of  magnesite.  Melts  were  deoxidized  by 
additions  of  0.5;  per  cent  aluminum.  Part  or  all  of  the  aluminum 
was  added  prior  to  tapping  at  3000°F.  In  all  cases  the 
temperature  of  the  melt  was  measured  using  an  optical  pyrometer 
prior  to  tapping.  No  further  temperature  measurements  were  made 
before  pouring,  j 

All  of  the  ingots  were  cast  in  plate  molds  1  inch  thick  by 
5  inches  wide  by  8  inches  high  topped  by  a  riser  which  tapered 
from  5  inches  by  1  inch  to  7  inches  by  3  inches  over  a  distance 
of  4  inches.  A  schematic  diagram  of  this  model  is  shown  in 
Figure  3-1.  The  mold  was  a  composite  of  CC>2  sand  and  exothermic 
molding  material.  The  exothermic  sleeve  was  1/2"  thick  at  the 
chill  and  had  a; 5/24  in/in  taper.  The  exothermic  material 
ignites  and  heats  to  a  temperature  above  that  of  the  melting 
point  of  steel.  Either  of  two  commmercial  exothermic  materials 
was  used:  Exomold  E  or  LD  Exothermic  Material.  To  avoid  gas 
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pick-up  the  cope  sections  were  all  baked  for  a  minimum  of  24 
hours  at  600°F.  The  drag  portion  of  the  mold  consisted  of  a 
water  cooled  copper  chill  block  which  was  given  a  light  zircon 
wash  and  supported  by  CO^  sand. 

Although  the  casting  weighed  20  pounds,  40  pound  heats 
were  always  melted  in  order  to  improve  chemical  and  temperature 
control.  The  iron-phosphorous  alloy  was  melted  under  argon  gas 
to  minimize  phosphorous  oxidation.  All  of  the  heats  were 
skimmed,  killed,  tapped  at  3000°F,  and  then  top  poured  directly 
into  the  mold.  The  exothermic  mold  was  not  ignited  before 
pouring.  Hot  water  was  first  run  through  the  chill  until  the 
molten  metal  was  poured  into  the  mold  at  which  time  cold  water 
was  used.  This  procedure  prevented  condensation  on  the  chill. 
After  filling  the  mold,  the  casting  was  capped  with  additional 
exothermic  material  to  reduce  top  heat  loss.  Structures  of 
the  castings  poured  were  fully  columnar,  with  the  columnar 
grains  extending  from  the  chill  up  into  the  riser.  Typical 
microstructures  are  shown  later  in  this  report. 

Thermal  measurements  were  made  during  the  solidification  and 
cooling  rate  of  one  of  the  iron-26  per  cent  nickel  ingots  cast, 
and  reported  previously1.  Twenty  mil  platinum-platinum-10  per 
cent  rhodium  thermocouples  were  inserted  horizontally  through 
the  exothermic  sleeve  in  the  5-inch  direction  at  distances  from 
the  chill  of  1  inch,  2  inches,  3  inches,  4  inches,  5  inches  and 
6  inches.  Each  thermocouple  was  protected  by  1/16-inch  I.D.  and 
1/8-inch  I.D.  fused  silica  tubes.  In  addition,  a  3/16-inch  fused 


silica  protection  tube  was  used  where  the  thermocouple 
penetrated  the  exothermic,  and  a  layer  of  CO^  sand  about  1/4- 
inch  thick  separated  each  thermocouple  ensemble  from  the 
exothermic  material.  Details  of  the  procedure  have  been 
described"*- . 

Cooling  curves  measured  for  the  unidirectionally  solidified 
ingot  are  shown  in  Figure  3-2.  The  sharp  change  in  slope  of 

i 

each  curve  indicates  the  time  the  dendrite  tips  reach  the  posi¬ 
tion  of  the  thermocouple  (due  to  the  higher  thermal  conductivity 
of  the  solid) .  It  is  clear  that  the  solidification  front 
progresses  unidirectionally  from  the  chill.  The  breaks  measured 
at  thermocouples  3,  4,  5  and  6  occur  at  very  close  to  the  same 
temperature,  namely  146 8 °C. 

From  these  cupves ,  the  position  of  any  isotherm  as  a  function 

of  time  is  readily  plotted.  From  such  plots,  the  position  of  the 

liquidus  and  solidus  were  shown  to  conform  closely  to  the 
1 

expressions  : 


.153 

/t7*  - 

.239 

(3 

”1) 

.145 

/P"  - 

.238 

(3- 

-2) 

where  t'  =  time  from  pour 

[ 

where  XL  is  position  of  liquidus  isotherm  (1468°C)  at  time  t', 

| 

and  Xg  is  position  of  solidus  isotherm  (liquidus  minus  25°C)*. 

Time  is  in  seconds  and  distance  in  inches. 

*  The  solidification  range  clearly  must  be  between  the  equilibrium 
range  (13°C)  and  the  range  for  no  solid  diffusion  (32°C). 
Reasoning  backwards  from  segregation  measurements  to  be  reported 


33 


The  "local  solidification  time",  t^,  as  a  function  of 
distance,  is  readily  calculated  from  equations  (3-1)  and  (3-2) 
as ; 

tf  =  4.84  X2  +  2.22  X  +  0.25  (3-3) 

and  averaging  cooling  rate,  z,  in  °C/sec  is  25 °C  divided  by  tf: 

e  =  (.194  X2  +  .089  X  +  .01)_1  (3-4) 

Cooling  rates  and  dendrite  arm  spacings  at  different  locations 
in  this  ingot  are  summarized  in  Table  4-1. 

B.  Other  Castings 

To  establish  relationships  between  dendrite  arm  spacing  and 
cooling  rate  (given  in  Chapter  4),  data  were  correlated  from 
earlier  studies  as  discussed  below. 

1 .  Furnace-cooled  iron-26%  nickel  alloy 

-3 

Data  for  this  alloy  in  the  cooling  rate  range  of  10  °C/sec 

to  5°C/sec  were  obtained  by  F.  R.  Mollard  and  M.  C.  Flemings  in 

earlier  research  supported  by  U.S.  Army  Materials  Research  Agency 
2 

at  M.I.T.  (during  the  period  December  9,  1961  to  December  8, 
1962) .  These  relatively  slow  cooling  rates  were  obtained  in  a 
controlled  solidification  furnace.  The  data  are  summarized  in 
Table  4-1. 

2 .  "Splat-cooled"  iron-26%  nickel  alloy 

In  a  separate  study  supported  by  United  States  Steel 
.  3 

Corporation  ,  small  droplets  of  iron-26%  nickel  alloy  were 
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levitation  melted  and  then  rapidly  solidified  ("splat-cooled") 
by  dropping  them  between  an  electrically  actuated  "hammer  and 
anvil".  Cooling  rates  were  measured  by  an  optical  sensor 
embedded  in  the  stationary  platen.  Results  are  summarized  in 
Table  4-1. 
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TABLE  3-1 


Heat  _ Nominal _  _ Actual* 


Number 

%  Ni 

% !  c 

%  Cu 

%  P 

%  Ni 

%  C 

%  Cu 

%  P 

%  A1 

1 

10.0 

i 

T 

_ 

. 

9.86 

_ 

0.06 

2 

15.0 

- 

- 

- 

14.7 

- 

- 

- 

0.04 

3 

20.0 

- 

- 

- 

19.9 

- 

- 

- 

0.21 

4 

26.0 

r 

- 

- 

26.0 

- 

- 

- 

0.24 

5 

26.0 

.12 

- 

- 

25.0 

0.12 

- 

- 

0.14 

6 

26.0 

.  33 

- 

- 

24.1 

0.33 

- 

- 

0.24 

7 

26.0 

.42 

- 

- 

24.6 

0.42 

- 

- 

0.18 

8 

26.0 

- 

- 

25.8 

- 

- 

- 

0.47 

9 

- 

■j- 

- 

4.0 

- 

- 

-4 

- 

10 

_ 

— 

25.0 

_ 

— 

— 

25.0 

- 

0.005 

*  Chemical  analyses  were  taken  at  approximately  1/2  inch  from 
the  chill.  0.5  per  cent  aluminum  added  as  deoxidizer  to  all 
heats. 


1500 


(Do)  3dfU\7H3dlAI3± 


Figure  3-2.  Cooling  curves  from  the  unidirectionally  solidified  ingot 
with  thermocouples  at  1.13,  2.17,  3.15,  4.17,  5.36,  and 
6.53  inches  from  the  chill. 
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CHAPTER  4  -  MORPHOLOGY  AND  DENDRITE  ARM  SPACINGS 

A.  Morphology 

Figures  4-1  to  4-6  show  selected  photomicrographs  from  this 
study ,  for  iron-nickel,  iron-nickel-carbon,  iron-phosphorous 
and  iron-copper  alloys. 

General  observations  from  these  photomicrographs  are  the 
following : 

(1)  The  dendrite  structures  tend  to  become  more  fully 
developed  with  increasing  distance  from  the  chill  (and 
hence  with  increasing  "local  solidification  time"  or 
"decreasing  cooling  rate").  Note,  for  example,  the 
horizontal  sections  of  the  iron-26  per  cent  nickel 
alloy  (Figures  4-1  and  4-2) . 

(2)  In  the  iron-nickel  alloys,  the  leaner  alloys  (10%  and 
15%  nickel)  have  less  fully  developed  dendrites  than 
the  richer  alloys.  In  the  transverse  sections  of 
Figure  4-3,  for  example,  the  lower  solute  containing 
alloys  appear  almost  "cellular".  We  continue  to  call 
these  dendrites,  however,  since  from  orientation  of  the 
transformation  structure,  it  is  clear  that  they  are 
growing  in  the  dendrite  direction,  <100>,  and  use  the 
term  "fibrous  dendrites"  when  no  side  branching  is 

visible , * 

*  We  prefer  to  distinguish  between  "cells"  and  "dendrites",  not 
by  presence  or  absence  of  secondary  dendrite  arms,  but  only  by 
the  criterion  that  the  dendrite  must  grow  in  or  near  its 
crystallographic  growth  direction. 
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(3)  Interstices  between  secondary  dendrite  arms  tend  to 
fill  in  as  previously  described  to  form  primary  plates 
(e.g.,  iron-26%  nickel)  or  continuous  "fibrous"  shapes 
(iron-10%  nickel,  iron-15%  nickel) , 

(4)  Iron-nickel  alloys  containing  carbon  are  more  readily 
etched  than  those  without  carbon  and,  perhaps  partly 
for  this  reason,  dendrites  in  the  ternary  alloys  appear 
better  developed  than  those  of  the  iron-nickel  binary 


(Figure  4-4) . 

(5)  In  iron-4  %  phosphorous  and  iron-25%  copper  alloys, 
horizontal  "secondary  plates"  develop  in  addition  to 


(the  less 


direct)  primary  plates.  Figures  4-5  and  4-6. 


B .  Dendrite  Arm  Spacing 


Qualitative  comparison  of  the  structures  of  Figures  4-1  to 

4-6  indicates  only  modest  effect  of  alloy  content  on  secondary 

dendrite  arm  spacing.  For  example,  for  a  location  2"  from  the 

chill,  secondary  dendrite  arm  spacings  in  all  alloys  studied 

were  between  75  and  90  microns.  Detailed  measurements  of  the 

1  2 

spacings  have  been  reported  '  . 

Alloy  analysis  influenced  primary  arm  spacing  much  more  than 
it  did  secondary  a;rm  spacing.  For  example,  in  Figure  4-4,  note 
the  primary  spacings  of  the  alloys  containing  .12%  carbon  and 
,42%  carbon  are  significantly  larger  than  those  of  the  two  alloys 
containing  0%  carbon  and  .33%  carbon.  Secondary  arm  spacings 
are,  however,  little  changed  by  the  change  in  analysis. 
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Dendrite  arm  spacing  is  strongly  influenced  by  distance 

from  the  chill  (and  hence  "local  solidification  time"  or 

"cooling  rate").  To  illustrate  this,  data  on  iron-26%  nickel 

3  4 

alloy  have  been  assembled  from  current  and  earlier  work  '  in 
Table  4-1  (as  described  in  Chapter  3).  These  same  data  are 
plotted  in  Figure  4-7  as  dendrite  arm  spacing  versus  cooling 
rate.  Note  that  over  a  range  of  cooling  rates  from  about 
10_3°C/sec  to  106°C/sec,  the  data  conform  approximately  to  a 
straight  line  (on  log-log  scales)  whose  slpe  is  about  one-third. 

For  comparison  with  the  structures  of  this  alloy  in  the 

unidirectional ly  solidified  plates.  Figures  4-8  and  4-9  show 

the  dendrite  structure  at  slower  and  faster  cooling  rates. 

Figure  4-8  shows  a  furnace  cooled  sample  (.44°C/sec  cooling 

rate)  and  Figure  4-9  shows  the  structure  of  two  splatted  samples 
S  6 

(at  2  x  10  and  2  x  10  °C/sec  cooling  rates) . 

C .  Discussion 

The  importance  of  "coarsening"  or  "ripening"  in  determining 
dendrite  arm  spacing  has  been  the  subject  of  a  number  of  recent 

rr  _  O 

studies  and  these  will  not  be  reviewed  here.  In  brief,  these 
studies  strongly  suggest  that  final  dendrite  arm  spacing  over  a 
very  wide  range  of  cooling  rates  is  determined  primarily  by  the 
rate  of  disappearance  of  dendrite  arms  during  solidification, 
not  by  the  number  that  initially  form.  Driving  force  for  this 
disappearance  of  dendrite  arms  is  reduction  of  liquid-solid 


surface  energy. 
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Calculations  on  this  ripening  (based  on  a  very  simple 
model)  indicate  the  critical  time  for  disappearance  of  a 
dendrite  arm,  i-s  proportional  to  the  following  group 

of  variables: 


HCl(1  -  k)  md' 


'crit 


(4-1) 


where  H  —  volumetric  heat  of  fusion 
CL  ~  liquid  composition 
k  =  equilibrium,  partition  ratio 
m  =  liquid  slope 

d  =  dendrite  arm  spacing  at  time  of  interest 
D  =  liquid  diffusion  coefficient 
T  =  temperature 


Insertion  of  reasonable  values  for  iron-26%  nickel  alloy 
in  equation  (4-1)  indicate  ripening  should  occur  in  these 
alloys  at  least  as  fast  as  in  other  alloys  for  which  it  has 
previously  been  shown  experimentally  to  occur.  The  linear 
dependency  in  Figure  4-7  of  dendrite  arm  spacing  on  cooling 
rate  is  strong  evidence  that  ripening  controls  the  final 
spacing  over  this  wide  range  of  cooling  rates . 


Note  that  it  should  not  be  surprising  that  ripening  can 
occur  at  the  extremely  short  times  involved  in  splat-cooling 
(or  chill  casting) .  Ripening  is  expected  to  be  significant 


when: 


(4-2) 


where  t^  =  local  solidification  time. 
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5-15 

Experimental  measurements  have  yielded,  for  a  variety 

of  alloys  (including  iron-26%  nickel.  Figure  4-7),  the 
relation: 

d^  =  Atf  =  A'/e  (4-3) 

where  df  =  final  dendrite  arm  spacing 
v , A , A '  -  constants 

e  =  cooling  rate 

Now,  from  (4-1)  and  (4-3)  ,  taking  d  as  a  constant  fraction  of  df: 

t  .  .  3-v 

«  d„  (4-4) 

tf  f 

Thus,  when  v  =  3,  a  constant  regardless  of 

cooling  rate,  and  ripening  can  be  expected  to  be  equally 
important  at  very  fast  cooling  rates  as  at  very  slow  cooling 
rates.  When  v  <  3,  extent  of  ripening  is  expected  to  be  more 
rapid  at  very  fast  cooling  rates  than  at  slower  ones. 

Most  experiments  on  cast  alloys  have  shown  v  to  be  between 
2  and  3,  and  in  these  cases,  therefore,  ripening  effects 
should  be  at  least  as  important  as  fast  as  at  slow  cooling  rates. 
Note  that  in  this  work  (Figure  4-7) ,  v  is  very  nearly  equal  to  3. 

Equation  (4-1)  is  certainly  an  oversimplified  picture  of  the 
ripening  process,  but  qualitative  conclusions  to  be  gained  from 
it,  such  as  the  above  considerations  of  effect  of  cooling  rate 
on  importance  of  ripening,  seem  justified.  Further,  the  general 
idea  that  final  observed  dendrite  structures  in  alloys  are 
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strongly  influenced  by  curvature  effects  (i.e.,  solidification) 
serves  to  explain  other  aspects  of  dendrite  structures.  As 
example,  the  tendency  of  dendrite  arms  to  fill  in  to  form 
"plates"  now  seems  clearly  to  result  from  this  effect. 

At  high  fraction  solids,  liquid-solid  surface  area  can  be 
reduced  more  efficiently  by  filling  in  spaces  between  rods  to 
form  plates  than  by  disappearance  of  dendrite  arms.  We  conclude 

that  when  plates  form,  this  is  the  primary  mechanism.  It  seems 

i 

| 

unlikely  that  heat  flow  alone  could  lead  to  these  plates  (as  has 

i 

been  previously  suggested)  since  they  are  often  observed  parti¬ 
cular  to  the  heat;  flow  direction  as  well  as  parallel  to  it 
(e.g.,  Figures  4-15,  4-6). 

D .  Conclusions  j 

(1)  Alloy  analysis  in  all  the  alloys  studied  had  relatively 
little  effect  on  secondary  dendrite  arm  spacing.  For  example, 
secondary  arm  spacings  2  inches  from  the  chill  in  unidirectionally 
solidified  ingots : were  between  75  and  90  microns  in  all  the 
alloys  (including  iron-niciel,  iron -nickel-carbon ,  iron- 

i 

phosphorous,  ironrcopper) . 

I 

(2)  In  iron-26%  nickel,  dendrite  arm  spacing  varies  linearly 
(on  log-log  scales)  with  cooling  rate  over  a  range  of  from 

_ 3  ^  i 

10  °C/sec  to  10  ^C/sec  cooling  rate.  Ripening  kinetics  control 
final  dendrite  arm  spacing  over  this  entire  range. 

(3)  Driving  force  for  formation  of  plate-like  structures  is, 
like  "ripening",  primarily  reduction  of  liquid-solid  surface  energy. 
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TABLE  4-1 


Experimental  Data  on  Cooling  Rates  and 
Dendrite  Arm  Spacings,  Iron-26  per  cent 
_ Nickel  Alloy _ 


Secondary  Dendrite 


Casting  Type 

Cooling  Rate 
°C/sec 

Arm  Spacing 
(microns ) 

Source 

Furnace-cooled 

4.62 

65 

(3) 

This  work  v 

II 

1.36 

X 

1 — 1 

1 

o 

p — 1 

77 

IV  II 

II 

6.35 

X 

I-* 

o 

J 

NJ 

145 

If  II 

If 

4.54 

X 

10"3 

280 

ii  II 

If 

4.36 

X 

ro 

1 

o 

1 — 1 

280 

It  Ii 

II 

4.28 

X 

10“3 

280 

II  VI 

H 

1.2 

X 

10-3 

605 

IV  VI 

Unidirectional 

Plate 

1/2"  from  chill 

9.6 

40 

This  Work 

2^  ii  ii  n 

3.4 

55 

It  M 

2  ii  »i  »• 

1.04 

80 

II  ll 

^  ii  H  ii 

2.9 

X 

i — 1 

i 

o 

1 — 1 

120 

IV  11 

Splat 

2.0 

X 

106 

.55 

(4) 

Strachan  v  1 

II 

2.0 

X 

10 1 6 

.60 

II 

II 

1.0 

X 

106 

.85 

II 

II 

5.0 

X 

105 

1.25 

IV 

II 

5.0 

X 

105 

1.00 

IV 

II 

3.3 

X 

105 

1.10 

ft 

n 

2.5 

X 

105 

1.8 

VI 

fi 

2.5 

X 

LO  1 

o 

1 — 1 

1.25 

11 

II 

2.0 

X 

iob 

1.25 

tv 

II 

1.0 

X 

105 

2.40 

11 

11 

8 . 0 

X 

104 

1.3 

VI 

11 

8.0 

X 

io4 

1.25 

11 
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Figure  4-1.  Iron-26  per  cent  nickel  alloy.  Photomicrographs  at  1", 
2",  and  4"  from  the  chill.  7.6X. 
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Figure  4-2.  Iron-26  per  cent  nickel  alloy.  Photomicrographs  at  1", 
2",  and  4"  from  the  chill.  34X. 
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Figure  4-3.  Iron-nickel  alloys  containing  10  per  cent  to  26  per  cent 
nickel.  Photomicrographs  at  2"  from  the  chill.  34X. 
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Figure  4-4.  Iron-26  per  cent  nickel  alloys  containing  0  per  cent 
carbon  to  0.42  per  cent  carbon.  Photomicrographs  at 
2"  from  the  chill.  34X. 
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Figure  4-5. 


Iron-4  per  cent  phosphorous  alloy, 
at  1",  2",  and  4"  from  chill.  34X. 


Photomicrographs 
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Figure  4-6.  Iron-25  per  cent  copper  alloy.  Photomicrographs  at  1", 
2",  and  4"  from  the  chill.  64X. 


FURNACE  COOLED 
UNIDIRECTIONALLY  SOLIDIFIED 
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Three-dimensional  view  of  iron-26  per  cent  nickel 
sample  furnace  cooled  at  4.3  x  10“-^°C/sec  . 


Figure  4-8. 
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Figure  4-9.  Two  splat  cooled  samples,  lron-26  per  cent  nickel  alloy^. 
Top-  2  x  10^°C/sec  cooling  rate,  1000X. 

Bottom-  2  x  105°C/sec  cooling  rate,  500X. 
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CHAPTER  5  -  MICROSEGREGATION 


A .  Procedure 

Chemical  segregation  existing  in  the  iron-nickel  and  iron- 
phosphorous  alloys  was  determined  by  electron  microprobe 
analysis.  Traces  were  made  on  samples  taken  in  a  plane  trans¬ 
verse  to  the  primary  growth  direction.  Maximum  and  minimum 
solute  contents  were  measured  and  "segregation  ratios" 
calculated . 

In  addition,  volume  per  cent  eutectic  in  the  iron-phosphorous 
system  was  measured  by  point  counting  in  the  following  way.  The 
polished  and  etched  structure  was  projected  on  the  ground  glass 
screen  of  the  metallograph  at  60X.  A  grid  of  63  points  was  placed 
over  the  field  of  view  and  the  points  that  fell  on  the  phase  were 
counted.  Between  30  and  60  grid  positions  were  counted  for  each 
sample  so  that  a  total  of  2000  to  4000  counts  were  amassed.  The 
ratio  of  points  falling  in  the  eutectic  was  compared  to  total 
number  of  points  yielded  per  volume  fraction  in  the  second  phase. 
Error  was  estimated  as  described  by  Hilliard  and  Cahn1  and  is 
shown  in  Table  5-7. 

B .  Iron-nickel  Binary  Alloys 

For  comparison  with  later  work,  Tables  5-1  and  5-2  summarize 

some  results  on  microsegregation  in  iron-26%  nickel  alloy 

obtained  in  a  previous  phase  of  this  research  and  reported 
2 

earlier  .  Small  samples  were  cooled  at  controlled  rate  in  an 
atmosphere  furnace.  Results  showed  the  following: 


1.  Segregation  ratio  is  not  strongly  affected  by  cooling 
rate  (Table  5-1)  . 

2.  Significant  diffusion  occurs  in  the  alloy  during  and 
after  solidification.  This  was  shown  by  rapid  quenching 
of  samples  from  at,  and  just  below  the  solidus  tempera¬ 
ture  (Table  5-2) .  Even  samples  quenched  from  just  at  or 
above  the  solidus  showed  much  less  segregation  than  would 

I 

be  predicted  by  the  "non-equilibrium  lever  rule". 

Table  5-3  summarizes  results  of  detailed  microprobe  studies  on 
unidirectionally  solidified  iron-26  %  nickel  and  iron-10%  nickel 

3 

alloys,  a  portion;  of  which  have  been  presented  earlier  .  Segrega¬ 
tion  ratio  decreases  slightly  with  increasing  distance  from  the 

i 

chill  in  iron-26%!  nickel,  but  not  at  all  in  iron-10%  nickel. 

I 

Solidification  times  at  these  various  locations  were  determined 
as  described  in  Chapter  3,  and  dendrite  arm  spacings  measured 
from  photomicrographs  such  as  those  in  Chapter  4 .  Results  are 
listed  in  Table  5^4. 

I 

I 

I 

Experimental  and  calculated  segregation  ratios  are  compared 
in  Table  5-5  for  these  two  alloys.  The  "calculated"  ratios  were 

I 

obtained  in  the  following  way: 

(1)  Values  of  ;ri  were  calculated  from  data  of  Table  5-4,  where 
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tf  -  solidification  time 
d  =  dendrite  arm  spacing  (primary  or  secondary) 
g  =  correction  factor  (as  given  in  footnotes  of  Table  5-4) 


(2)  Calculated  segregation  ratio  at  room  temperature ,  S°, 
were  then  determined  directly  from  the  appropriate  figure  of 

Chapter  2 . 


The  correction  factors,  g,  applied  to  the  dendrite  arm 
spacing  were  arbitrary  except  that  the  same  correction  factor 
was  used  for  both  alloys  and  all  solidification  times.  The 
factor  chosen  was  that  which,  for  a  given  model,  gave  closest 
agreement  with  experiment.  These  were: 


Plate  model,  primary  spacing 
Plate  model,  secondary  spacing 
Cylinder  model,  primary  spacing 
Cylinder  model,  secondary  spacing 


g  =  0.13 

g  =  0.30 

g  =  0.17 

g  -  0.40 


A  correction  factor,  g,  of  1  would  indicate  direct 
correspondence  of  experiment  with  theory.  Factors  less  than  1 
indicate  considerably  less  segregation  than  that  predicted.  Note 
that  the  cylindrical  models  give  closer  correspondence  with 
experiment,  but  the  required  correction  factor  is  still  signifi¬ 
cantly  less  than  1.0.  This  is  discussed  further  below. 


C .  Iron-nickel-carbon  Alloys 

Carbon  was  added  to  an  iron-2 6%  nickel  alloy  in  amounts  up 
to  0.42%,  and  segregation  ratios  measured.  Table  5-6.  No 
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significant  effect  was  measured. .  A  result  considerably 

I 

different  from  that  is  obtained  when  carbon  is  added  to  iron- 
chromium  alloy^. 


D .  Iron-4 %  Phosphorous  Alloy 

Results  of  lineal  analysis  of  the  iron-4%  phosphorous  alloy 
are  listed  in  Table  5-7.  The  per  cent  eutectic  decreases  only 
slightly  with  increasing  distance  from  the  chill.  Moreover, 
using  solidification  times  and  dendrite  arm  spacings  listed  in 
Table  5-8,  results  agree  with  experiment  assuming  correction 
factors,  g,  are  not  greatly  different  than  those  for  the  iron- 
nickel  alloys  (g  -  0.19  for  primary  spacing,  plate  model; 
g  =  .56  for  secondary  spacing,  plate  model).  Results  are  listed 
in  Table  5-9. 

Table  5-10  lists  minimum  solute  contents  measured  with  the 
microprobe.  In  agreement  with  experiments  on  other  alloys, 
these  minima  are  also  relatively  constant  with  increasing  distance 
from  the  chill.  However,  they  are  substantially  higher  than  would 
be  predicted  from  Figure  2-11  using  correction  factors,  g,  such 
as  have  been  used  above.  Much  smaller  factors  (g  <  0.10)  must  be 
employed.  The  reason  for  the  higher  solute  content  than  expected 
(and  hence  lower  required  g)  is  not  known. 


E.  Discussion 


When  early  work  was  conducted  on  microsegregation  in  ferrous 
and  non-ferrous  alloys,  it  was  observed  that  correction  factors 


of  about  0,2  to  0 


3  (applied  to  dendrite  arm  spacing)  were 
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necessary  to  make  experiment  and  theory  correspond.  At  that 
time,  only  a  single  plate  model  was  employed,  and  the  procedure 
employed  was  valid  only  for  relatively  small  amounts  of  diffu¬ 
sion  in  the  solid  during  and  after  solidification. 

Work  reported  in  Chapter  2  of  this  report  was  therefore 
undertaken  to  permit  better  correspondence  of  theory  with 
experiment.  A  computer  analysis  valid  for  extensive  diffusion 
was  developed.  Also,  a  cylindrical  dendrite  model,  as  well  as 
plate-like  dendrite  model,  was  examined.  Agreement  with  experi¬ 
ment,  as  noted  above,  was  better.  For  iron-nickel  alloys,  the 
correction  factor  rose  as  high  as  0.4  (for  cylindrical  model, 
secondary  spacing)  but  this  factor  is  still  too  far  from  1  to 
be  explained  by  considerations  such  as  inaccurate  diffusion 
data,  more  complex  dendrite  geometry,  etc. 

The  discrepancy  is  resolved  only  in  the  light  of  recent 
experiments,  mentioned  previously  in  Chapter  4,  which  show  that 
substantial  "ripening"  or  "coarsening"  occurs  during  dendritic 

5-9 

solidif ication  .  In  essence,  the  neglect  of  curvature  effects 
is  now  shown  to  be  an  imperfect  approximation.  In  fact,  the 
curvature  effects  mean  that  throughout  dendritic  solidif ication , 
small  arms  are  disappearing,  and  larger  ones  growing  at  their 
expense.  The  dendrite  arm  spacing  finally  observed  in  a 
casting  or  ingot  is  not  that  existing  throughout  freezing  but  is 
much  coarser  than  that  that  was  present  relatively  early  in 
freezing . 
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Thus,  the  correction  factor,  g,  reflects  in  part  at  least, 
the  fact  that  average  dendrite  arm  spacing  during  solidification 
is  significantly  smaller  than  that  finally  observed.  Considerable 
solid  diffusion,  previously  unaccounted  for,  is  taking  place  in 
small  dendrite  arms  which  are  "remelting"  and  re-precipitating  on 
the  larger  arms  throughout  solidification. 

F.  Conclusions 

1.  Segregation  ratio  in  iron-26%  nickel  alloy  decreases  only 
slightly  with  increasing  distance  from  the  chill.  Segregation 
ratio  in  iron-10%  nickel  alloy  is  essentially  constant. 

I 

2.  In  iron-4%  phosphorous  alloys,  amount  of  second  phase  and 

I 

minimum  solute  content  are  essentially  constant  with  increasing 
distance  from  the  chill. 

3.  Carbon  in  amounts  up  to  0.42%  does  not  significantly 
affect  segregation  ratio  of  iron-26%  nickel  alloy. 

4.  Good  agreement  of  calculated  segregation  with  that 
experimentally  observed  in  iron-nickel  alloys  is  obtained, 
provided  calculations  are  based  on  measured  dendrite  arm  spacing 
times  a  correction1  factor,  g.  The  required  correction  factor 
depends  on  dendrite  model  employed  but  not  on  solidification  time 
or  alloy  content.  j Depending  on  the  model  employed,  g  lies  between 
0.17  and  0.40.  i 
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5.  It  appears  that  the  need  for  a  correction  factor,  g, 
arises  primarily  because  of  the  "coarsening"  or  "ripening"  now 
known  to  occur  during  dendritic  solidification. 

6.  In  agreement  with  theory,  segregation  ratio  of  iron-10% 
nickel  alloy  is  greater  than  that  in  iron-26%  nickel  alloy. 
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Table  5-1 

Furnace  Cooled  Iron-25%  Nickel  Alloy* 


Dendrite  Arm 


Cooling  Rate 
( °C/sec) 

Spacing 

(d) 

CM 

cra 

S° 

.14 

77 

28.2 

22.8 

1.23 

-9 

6.4  x  10 

145 

28.7 

22.8 

1.25 

-3 

4.3  x  10 

280 

30.7 

26.7 

1.15 

* 


Samples  were  water-quenched  55°C  below  the  liquidus  temperature. 


Table  5-2 


6  6 

Segregation  Ratio  in  Interrupted  Cooling  Experiments 
(Cooling  Rate  4.3  x  10  ^  °C/sec) 


Quenching 

Temperature 

CM 

C 

m 

VCm 

Tm  -  21 °C 

30.2 

23.5 

1.29 

tm  -  25°C 

i 

29.7 

25.5 

1.17 

T  -  55 °C 

M 

30.7 

26.7 

1.15 

Note:  Non-equilibrium  solidus  is  approximately  - 

(see  footnote  to  equation  3-2). 


25  °C 
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Table  5-3 


Electron  Probe 

Microanalyzer 

Results 

S° 

e 

c°/c° 

M_ _ rn. 

Alloy  -  Distance 

from  chill 
( inches) 

C° 

M 

(%  Ni) 

C° 

m 

(%  Ni) 

Fe-26%  Ni 

ii  ILL 

1/2 

32.6 

25.5 

1.28 

1 

30.0 

24.5 

1.26 

2 

29.5 

25.2 

1.17 

4 

29.8 

25.2 

1.18 

Fe-10%  Ni 


1/2 

11.6 

8.59 

1.35 

1 

11.3 

8.58 

1.32 

2 

11.6 

8.40 

1.38 

4 

11.6 

8.40 

1.38 
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Table  5-4 

Measured  Dendrite  Arm  Spacing  and  Solidification  Times 


in  Iron 

-Nickel  Alloys 

Alloy 

^Distance 
from  chill 
(inches ) 

Primary 
dendrite  arm 
spacing,  d 
(microns) p 

Secondary 
dendrite  arm 
spacing,  d 
(microns ) 

Solidi¬ 
fication 
time,  tf 
( seconds ) 

Fe-26% 

Ni 

|  1/2 

70 

40 

2.6 

Fe-26% 

Ni 

1 

100 

55 

7.5 

Fe-26% 

Ni 

2 

190 

80 

24 

Fe-26% 

Ni 

4 

360 

120 

86 

Fe-10% 

Ni 

!  1/2 

50 

- 

2 . 6 

Fe-10% 

Ni 

1 

;  i 

75 

- 

7.5 

Fe-10% 

Ni 

2 

125 

80 

24 

Fe-10% 

Ni 

,  4 

235 

150 

86 
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Table  5-5 

Comparison  of  Measured  and  Calculated  Segregation 
Ratios  in  Iron-Nickel  Alloys 

Exp.  Calculated  Segregation 


Alloy 

Location 
(inches 
from  chill) 

Measured 
Segregation 
Ratio,  S° 

Ratios , 

S° _ 

(1) 

(2) 

(3) 

(4) 

Fe-26%  Ni 

1/2 

c 

1.28 

1.29 

1.42 

1.29 

1.32 

Fe-26%  Ni 

1 

1.22 

1.22 

1.21 

1,21 

1.29 

Fe~2 6 %  Ni 

2 

1.17 

1.17 

1.23 

1.23 

1.22 

Fe-26%  Ni 

4 

1.18 

1.18 

1.23 

1.23 

1.12 

Fe-10%  Ni 

1/2 

1.35 

1.35 

- 

1.58 

- 

Fe-10%  Ni 

1 

1.32 

1.32 

- 

1.35 

- 

Fe-10%  Ni 

2 

1.38 

1.38 

- 

1.33 

- 

Fe-10%  Ni 

4 

1.38 

1.38 

- 

1.33 

- 

(1)  Based  on  plate  model,  using  primary  dendrite  arm  spacing  times 
0.13  to  calculate  ri . 

(2)  Based  on  plate  model,  using  secondary  dendrite  arm  spacing 
times  0.30  to  calculate  n* 

(3)  Based  on  cylinder  model,  using  primary  dendrite  arm  spacing 
times  0.17  to  calculate  n. 

(4)  Based  on  cylinder  model,  using  secondary  arm  spacing  times 


0.4  to  calculate  n . 
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Table  5-6 

The  Effect  of  Carbon  on'  Microsegregation 
in  an  Iron-26%  Nickel  Alloy  -  2  Inches 
_ from  the  Chill 


Chemical 

Analysis  C°„ 

M 

Wt.  %  C  1  Wt.  %  Ni 


C° 

m 

Wt.  %  Ni 


S° 

e 


0 


29.5  25.5  1.17 


.33 


26.8  22.4  1.19 


42 


26.4 


22.7 


1.16 
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Table  5-7 

Volume  Per  Cent  Eutectic  Measured  Metallographically 
at  Different  Distances  from  the  Chill 


Distance  from 
the  Chill 


1" 


2" 


3  " 


4" 


Volume  Per  Cent 
28.3  +  1.5* 

28.5  +  1.5* 

27.5  +  1.5* 


Weight  Per  Cent 

26.4 

26.5 

25.7 


27.0  +  1.5* 


25.2 
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Table  5-8 

Measured  Dendrite  Arm  Spacings  and  Solidification  Time 
_ in  Iron-4  %  Phosphorous  Alloy _ 


Distance 
from 
chill 
( inches ) 


1 


Primary 
,  dendrite* 
arm  spacing, 
_d  ,  (microns) 


122 


Secondary 

dendrite  Solidification 

arm  spacing ,  time ,  tf 

d  ,  (microns)  (seconds) 


50 


7.5 


220 


73 


24 


300 


100 


51 


380 


108 


86 


73 


Table  5-9 

Comparison  of  Measured  and  Calculated  Weight  Per  Cent 
Eutectic  in  Iron-4, 0  Per  Cent  Phosphorous  Alloy 


Location 

(inches 

from 

drill) 


Experimentally 
Measured  Weight 
Per  Cent  Eutectic 


Calculated  Weight 
Per  Cent  Eutectic 


(1)  Based  on  plate  model,  using  primary  dendrite  arm  spacing 
times  0.19  to  calculate  r) . 

(2)  Based  on  plate  model  using  primary  dendrite  arm  spacing 
times  0.56  to  calculate  n. 
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Table  5-10 

Electron  Microprobe  Analysis  for  Minimum  Phosphorus 
Contents  at  Different  Distances  from  the  Chill 


Distance  from  Minimum  Solute  Content 

the  chill  C  (wt.  %  P) 

X ,  (inches) _  _  m 


1/2  |  2.2 

2  |  2.1 


4 


2.0 


APPENDIX  A:  List  of  Symbols 


Composition  of  the  liquid  within  the  volume 
element 

Local  maximum  solid  composition  at  non-equilibrium 
solidus  and  at  room  temperature,  respectively 

Local  minimum  solid  composition  at  non-equilibrium 
solidus  and  at  room  temperature,  respectively 

Solid  composition  within  the  dendrite 

Solid  composition  in  the  j  slice  and  at  the  k 
interval  of  time 

Solid  composition  at  the  liquid  solid  interface 
Dendrite  arm  spacing  at  time  of  interest 
Final  dendrite  arm  spacing 

Diffusion  coefficient  in  liquid  and  solid  phases, 
respectively 

Weight  fraction  eutectic,  liquid,  solid, 
respectively 

Microsegregation  correction  factor 
Volumetric  heat  of  fusion 
Partition  ratio 

One  half  the  dendrite  arm  spacing 
AA2/DsAt 

Slope  of  the  liquidus 

Segregation  ratio,  C M/C  ,  at  the  non-equilibrium 
solidus  and  at  room  ^temperature ,  respectively 

Time  from  the  beginning  of  solidification  within 
the  volume  element 

Local  solidification  time 

Time  from  beginning  of  solidification  of  the  ingot 
Critical  time  for  disappearance  of  a  dendrite  arm 
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T 

X 

V  Xs 

a 

e 

X 

A. 

1 

n 


Absolute  temperature 
Distance  from  the  chill 

Position  of  liquidus  and  solidus,  respectively 

Dstf/L2 

Cooling  rate 

Distance  from  the  dendrite  centerline 

Position  of  the  solid  liquid  interface  within  the 
volume  element 

tf/L2 
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APPENDIX  B:  Method  III  -  Finite  Difference. 

Mass  Balance  Technique 

The  major  assumptions  of  the  solidification  model  are 
presented  at  the  beginning  of  Chapter  2 .  First  consider  the 
dendrites  to  be  plate-like  and  the  characteristic  volume  element 
to  be  as  in  Figure  B-l.  The  initial  condition  for  the  computa¬ 
tion  is 


when  t  =  0,  A.  =  0 ,  C  ( A  •  )  -  kC  ,  C  =  C  (B-l) 

X  o  X  V  ±j  V 

The  mechanics  of  the  numerical  technique  are  illustrated  in 
Figure  B-l  and  are  described  below. 


The  Solidification  Interval 

(1)  The  composition  distribution  in  the  solid  of  a  plate-like 
dendrite  is  schematically  represented  by  the  heavy  line  (C  ) 
in  Figure  B-l  where  the  subscript  j  indexes  the  position  of  the 
slice  and  the  subscript  k  refers  to  intervals  of  time.  The 


composition  of  the  liquid  is  calculated  from  the  mass  balance 

-i  * 

C  A  A 


1  3 

c  -  —  y 
o  L  "  s  . 

D-l  J 


(B-2 ) 


where 


(B-3 ) 


and  j*  refers  to  the  slice  at  the  liquid  solid  interface.  (Note: 
For  the  first  step  C  =  C  ) 

Xj  q 
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(2)  The  interface  is  advanced  a  distance  A  A  by 

transformation  of  an  increment  of  solid  of  composition  C  ( A) 

Sk+1 

which  is  in  equilibrium  with  the  liquid  CT  as  given  by  the 

Lk 

phase  diagram. 


(3)  The  temperature  of  the  volume  element  is  determined 

from  the  phase  diagram  and  the  value  of  CT  . 

;  Lk 

(4)  Diffusion  is  allowed  to  occur  in  the  solid  for  an  amount 

of  time  At  which  is  the  time  elapsed  as  the  interface  moves  ahead 

one  slice  from  A,  to  A. 

I  Xk  ^+1 

Solid  diffusion  is  computed  by  iteratively  evaluating  the 
finite  difference  algarithm  for  Fick's  Second  Law 


C  +  (M  -  2)  C  +  C 

j+l,k-l _ S j  t k^l  5j-l,k-l 

|  M 

|  ..  (B-4) 


where 

M 

for  each  slice  within  the  solid,  1  £  j  £  j*.  The  elapsed  time, 
t,  may  be  computed  from  the  linear  or  parabolic  growth  expres¬ 
sions  (equations,  2-1  or  2-2)  or  from  thermal  data,  in  the  form 
of  an  expression  or  table,  that  relates  temperature  and  time. 

If  the  diffusion  coefficient,  D  ,  is  a  function  of  both 

!  S 


A  A 
DAt 


(B-5) 


*  H.  S.  Carslon,  and  J.  C.  Jaeger:  Conduction  of  Heat  in  Solids, 
Clarendon  Press,  Oxford,  1959,  p.  466. 
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temperature  and  composition,  it  must  be  re-evaluated  for  every 
iteraction  of  the  algarithm. 

The  finite  difference  approximation  is  most  accurate  for 

large  values  of  M  and  the  solution  will  be  unstable  for  values 

2 

of  M  less  than  two.*  For  this  reason  whenever  AA  /D  At  was 

less  than  4,  diffusion  was  considered  to  take  place  over 
2 

p  >  4AA  /D  At  time  increments  each  A9  in  length,  such  that 

2 

A  A  /D  A9  >  4  and  pA0  =  t.  The  solution  for  the  incremental 
time  period  At  is  then  attained  by  repeating  equation  C-4  over 
all  slices,  j,  p  times. 

The  boundary  conditions  for  this  computation  are 
3  C 

at  X  -  0 ,  — ■  =  0;  and  at  A  =  A.  ,  C  =  C  (A.)  (B-6) 

3  A  is  sk  i 

where  C  (A.)  is  computed  by  mass  balances  as  in  steps  (1)  and 

S  n  1 

k 

(2).  And  the  composition  of  the  first  slice  in  the  -A  direction 
is  always  set  equal  to  the  composition  computed  for  the  first 
slice  in  the  +A  direction.  Thus,  the  flux  across  A  =  0  is  zero* 

(5)  After  diffusion  has  occurred  for  a  time  At  and  a  new 
composition  distribution  has  been  computed  steps  (1)  through  (4) 
are  repeated;  a  new  liquid  composition  is  computed,  the  composi¬ 
tion  of  a  new  solid  increment  and  the  temperature  are  found  from 
the  phase  diagram  and  diffusion  is  again  allowed  to  occur.  The 

*  H.  S.  Carslon,  and  J.  C.  Jaeger:  Conduction  of  Heat  in  Solids, 
Clarendon  Press,  Oxford,  1959,  p.  466. 


above  procedure  is  repeated  until  the  alloy  is  all  solid  or 
until  a  eutectic  composition  is  reached  in  the  liquid. 


(6a)  Eutectic  Systems:  During  solidification  of  an  alloy 

i 

in  a  eutectic  system,  if  the  liquid  composition  reaches  the 
eutectic  composition  at  a  time  t^,  diffusion  is  computed  in  the 

primary  phase  at  the  eutectic  temperature  until  solidification 

i 

is  completed,  i.e..  At  =  t^  -  t£.  As  a  result  of  diffusion  of 
solute  into  the  primary  phase,  some  eutectic  must  disappear; 
and  this  shift  in  the  primary  phase-eutectic  interface  is 

j 

monitored  throughout  the  calculation.  The  following  boundary 
conditions  are  used  in  the  computation  from  time  tE  to  tf: 


A  =0,  ^  =0;  A  =  Ag ,  C*  =  Cg  (maximum  solubility)  (B-7) 


where  AE  represents  the  distance  to  the  eutectic  interface. 

l 

In  the  special  case  that  the  last  increment  of  solidifi¬ 
cation  is  reached  (i.e.,  A.  =  L  -  AA)  and  the  liquid 

i  i 

composition  is  less  than  the  eutectic  composition  and  greater 
than  the  maximum;  solubility  limit  some  eutectic  is  considered 
to  have  formed  and  the  amount  is  given  by 


f 


E 


CL  -  CQ  (maximum  solubility) 
CE  -  Cg  (maximum  solubility) 
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where  f„  =  weiqht  fraction  eutectic 
E 

C  =  composition  of  the  eutectic  at  the  eutectic 
E  temperature 

C  (maximum  solubility)  =  composition  of  the  primary 
&  phase  at  the  eutectic  temperature 

In  the  above  two  cases  the  non-equilibrium  solidus  temperature 
is  considered  to  be  the  eutectic  temperature. 

In  the  case  that  the  last  increment  of  solidification  is 
reached  and  the  liquid  composition  is  less  than  the  maximum 
solubility  limit,  the  computation  is  handled  in  the  same  manner 
as  a  single  phase  alloy. 

(6b)  Single  Phase  Alloy:  Since  the  compositions  gradients 
near  the  end  of  solidification  are  usually  steep,  the  average 
composition  of  the  last  increment  to  freeze  may  not  be  a  good 
indication  of  the  non-equilibrium  solidus  temperature.  It  is 
often  desirable  to  look  at  the  last  slice  to  freeze  in  more 
detail  (and  in  the  time  sharing  version  of  this  program  the 
decision  is  left  to  the  human  operator,  see  Appendix  D) .  When 
desirable,  the  last  increment  to  freeze  is  divided  into  smaller 
slices  and  the  composition  distribution  within  that  increment 
calculated  in  greater  detail.  The  last  increment  is  broken 
into  one  fifth  the  number  of  slices  that  were  used  in  the 
original  computation  and  the  analysis  repeated  for  the  new  C  , 

L f  and  0  .  The  factor  of  0.2  was  found  to  be  adequate  to 
describe  the  final  solidification  and  the  keep  computer  time 
reasonable.  The  composition  of  the  last  solid  to  freeze  in  this 
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final  calculation  is  then  used  to  define  the  non^equilibrium 
solidus  temperature. 

Cooling  from  the  Non-Equilibrium  Solidus 

Following  tjhe  completion  of  solidification,  the  change  in 
solute  redistribution  that  occurs  on  cooling  from  the  non¬ 
equilibrium  solidus  to  room  temperature  may  be  computed  using 
the  finite  difference  technique.  The  cooling  rate  relation  may 
be  input  to  the  computer  from  hypothetical  or  empirical  tabular 
data  or  in  the  form  of  an  expression.  The  cooling  rate  relation 
defines  the  time ; for  diffusion  within  each  temperature  interval 
and  the  diffusion  coefficient  is  calculated  on  the  basis  of  the 
average  temperature  within  the  temperature  interval.  The  amount 
of  diffusion  is  gomputed  by  evaluation  of  the  algorithm  (B-4) 
at  each  slice  j  within  each  temperature  interval  and  the  size  of 
the  temperature  interval  is  chosen  to  keep  M  greater  than  4.  The 
initial  condition  is  the  solute  distribution  at  the  solidus 
temperature.  The  boundary  conditions  for  a  single  phase  alloy 
are 

A  =  0,  =  0;  A  =  L,  |£  =  0  (B-9) 

and  for  a  eutectic  alloy 

X  =  0 ,  “S  =0;  A  =*  £  (1  -  f  )  ,  C*  =  C  (max.  solubility) 

o  A|  £j  S  S 

•  .  (B-10 ) 

Computations  are  continued  until  room  temperature  is  reached  or 
until  the  effect  of  diffusion  on  microsegregation  becomes 
negligible . 
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Cooling  Curves:  Since  diffusion  can  be  described  in 
terms  of  values  it  is  convenient  to  choose  a  one  second 
solidification  time  and  to  vary  the  dendrite  arm  spacing  in 
order  to  get  the  microsegregation  parameters  in  terms  of  n . 

If  actual  cooling  curves  are  to  be  used,  this  procedure 
necessitates  the  readjustment  of  the  data  to  produce  a  meaning¬ 
ful  result, 

Iron-Nickel :  To  evaluate  the  cooling  rates  for  various 

value  of  n,  the  six  thermocouple  measurements  taken  at  various 
distances  from  the  chill  in  an  iron-26%  nickel  alloy  and  des¬ 
cribed  in  Chapter  3  were  plotted  in  a  special  way.  First,  the 
solidification  range  was  taken  to  be  25°C.  The  time  to  drop 
25 °C  was  measured  from  each  cooling  curve  and  the  six  values 
defined  as  the  solidification  times,  t^.  Temperature  readings 
were  then  taken  from  each  curve  at  multiples  of  its  own  t^-  and 
plotted.  The  resulting  curves  overlapped  and  deviated  less 
than  +  15 °C  down  to  1150°C  and  approximately  18  Three 

positions  overlapped  +  10°C  all  the  way  down  to  1000°C.  This 
allowed  the  master  curve  of  temperature  versus  O^.to  be  made 
which  ranges  from  1469.5  to  1 000°C  and  is  shown  in  Figure  B-2. 
These  values  were  then  used  for  all  g's  in  the  iron-nickel 
system  as  follows: 


84 


(1)  The  temperature  of  the  zero  time  point  of  the  master 
curve  and  the  liquidus  point  for  the  alloy  were  made  to  coincide 
by  shifting  the: master  curve. 

(2)  The  solidification  temperature  range  for  a  given  n  is 
taken  from  the  computer  solution  after  solidification  has  occurred 
and  a  solidification  time  is  obtained  from  the  liquidus  and 
solidus  temperature,  and  the  cooling  curve.  The  solidification 
time  is  then  compared  with  the  assigned  value,  i.e,  1  second 

and  the  remaining  times  in  the  solidification  curve  adjusted  by 
this  factor.  ; 


Iron-Phosphorus :  In  the  iron-phosphorus  system  cooling 

rates  of  200°C/sec  and  800°C/sec  were  used  in  conjunction  with 
a  1  second  solidification  time. 


Cylindrical 


Geometry :  In  some  cases  a  cylindrical  volume 


element  may  be  of  more  interest  than  a  plate-like  volume  element. 
In  such  a  case  the  axis  of  the  cylinder  is  considered  to  coincide 
with  the  dendrite  spine  and  the  liquid  solid  interface  advances 
from  the  axis  at  jtime  t  *  0  to  a  radius  of  L  at  time  t^,  where 

i 

L  is  one  half  the  dendrite  arm  spacing  and  t^  is  the  solidifica¬ 
tion  time.  The  rate  of  advance  of  the  interface  may  be  given  by 
a  linear  volume  growth  rate  or  a  parabolic  volume  growth  rate: 


dV 

dt 

dV 

dt 


constant 

I 

constant 
✓t  1 


r . 
l 


r . 
i 


=  constant  /t 


(B-ll) 


=  constant  4/t 


(B-12 ) 
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where  denotes  the  position  along  the  radius  of  the  solid- 
liquid  interface.  In  addition,  as  with  the  plate-like  model, 
the  position  of  the  interface  may  be  calculated  from  a  cooling 
curve.  In  cylindrical  coordinates  Fick's  Second  Law  is  written 


dC 

3t 


D  ( 


92C  1  3C 
3r2  r  9r 


( B  — 1 3 ) 


which  is  transformed  by  the  finite  difference  method 


C  .  , 
5/k 


(1  +  i^-)  Cj  +  1  +  (M  -  2 )  C_.  +  (1 

M 


0.5 


)  C. 


3-1  J-1 


(B-14) 


The  mechanics  of  the  solution  are  the  same  as  that  used  for 
the  plate-like  model  except,  now,  equal  slices  are  taken  along 
the  radius,  the  algorithm  (B-14)  is  used  to  evaluate  the  diffu¬ 
sion  in  the  solid  and  the  rate  of  advance  of  the  interface  is 
given  by  expressions  such  as  (B-ll)  or  (B-12) . 


S(  J,k  +  1) 


Figure  B-l.  Schematic  diagram  of  the  sequence  of  steps 

in  the  Mass  Balance  Technique  for  calculating 
solute  redistribution  in  an  alloy  solidifying 
into  a  single  solid  phase  of  variable  com- 
i  position. 
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